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Deconvolution  methods  have  been  applied  to 
measurements  made  with  different  electrical  sensors 
including  resistive  and  capacitive  dividers. 
Deconvolved  and  directly  measured  waveforms  have  been 
compared  with  good  results. 

Introduction 

In  order  to  characterize  a  sensor,  we  use  linear 
system  methods  to  compute  the  system  step  (or  impulse) 
response.  This  involves  unfolding  a  convolution 
integral  equation  (see  (1)  below)  in  the  presence  of 
noise.  We  construct  the  response  by  digital 
computation.  Using  the  derived  response  and  a 
measured  output  waveform,  the  same  software  allows  us 
to  determine  the  associated  input  waveform. 
Comparison  of  the  calculated  input  with  the  measured 
input  shows  that  features  of  interest  can  be 
reconstructed  faithfully.  However,  no  systematic 
study  of  errors  is  presented. 

Description  of  the  Deconvolution  Method 

In  our  laboratory  an  input  signal,  x(t),  to  a 
device  and  the  corresponding  output  signal,  y(t),  from 
it  are  recorded.  For  a  system  which  is  linear, 
causal,  and  time-shift  invariant,  linear  systems 
theory  implies 

y(t)  =  x(s)  •  g(t-s)ds,  (1) 

where  g  is  the  impulse  response  of  the  system  [1]. 
Using  (1)  we  approximate  g  (or  its  integral,  the  step 
response,  h)  given  estimates  of  y  and  x. 
Alternatively,  we  determine  x  given  y  and  g  (or  h). 
Both  mathematically  and  computationally,  these  two 
problems  may  be  solved  by  identical  means  because  x 
and  g  can  be  interchanged  in  (1).  The  principal 
difficulty  in  solution  is  that  (1)  is  an  integral 
equation  of  the  first  kind  and  is  singular  or  nearly 
so. 

In  our  case,  x  and  y  are  observed  over  a  time- 
window,  OStSW,  and  consist  , of  discrete,  digitized 
measurements  which  contain  digitization  errors  and 
system  noise.  Thus,  each  of  the  vectors  x  =  (X|^), 

y  =  ®  ^  ^^k^  ®  time  series.  The 

subscript  k  =  1,  ...,  N  and,  with  the  measurement 
system  used,  N  must  be  ^  512,  When  discretized,  (1) 
gives  a  system  of  linear  equations  which  in  matrix 
form  is 

Xg  =  y,  (2) 

and  X  is  an  N  X  N  Toeplitz  matrix  with 


0  for  1ii<j5N 


Equation  (2)  may  be  inconsistent  and  have  no 
solution.  It  will  certainly  be  ill-conditioned  so  g 
is  sensitive  to  small  changes  in  y.  Precisely  because 
y  is  noisy,  g  is  not  accurately  determined  on  the 
basis  of  (2)  above.  In  particular,  solutions  tend  to 
be  highly  oscillatory  and  to  take  on  large  values; 
neither  is  physically  correct.  In  order  to  restore 
physical  meaning  to  the  mathematical  computation  of  g, 
we  change  the  problem.  Equation  (2)  is  altered  so 
that  simultaneously  neither  inconsistency  nor  ill- 
conditioning  dominates  the  solution. 

Instead  we  consider  the  related  minimization  [2]: 

rain  ( I  |Xg-yl  +  y|  |Cgl  f)  .  (3) 

II  •  II  denotes  the  root  mean  square  (rms)  norm,  7  is 
the  regularization  parameter,  and  C  is  a 
regularization  operator.  We  take  it  to  be  the  second 
symmetric  difference  operator  [3].  ||  Cg  ||  is  a 

measure  of  the  oscillation  in  g.  For  7=0,  the 
solution  to  (3)  agrees  with  that  to  (2)  if  one 
exists.  Thus,  in  this  case,  g  is  consistent  with  the 
data  X  and  y.  However,  as  we  have  seen  above,  mere 
data  consistency  typically  produces  estimates  of  g 
that  are  useless. 

For  7  >  0,  the  solution  to  (3)  is  no  longer 
consistent  with  the  measurements.  The  particular  form 
of  the  second  term  is  not  dictated  bythe  need  to 
stablize  the  computation  of  g.  Instead,  it  arises  in 
trying  to  incorporate  qualitative  information  that  is 
entirely  independent  of  the  physical  measurements.  In 
our  case,  this  prior  information  is  based  on 
experience  with  the  characteristics  of  the  sensors 
under  study;  the  circuits  involved  have  a  frequency 
response  which  does  not  allow  for  a  highly  oscillatory 
step  response.  By  minimizing  the  sum  of  these  two 
terras,  we  find  a  solution  which  simultaneously  comes 
close  to  being  data  consistent  and  is  not  physically 
unreasonable. 

The  quadratic  minimization  (3)  is  equivalent  to 
the  linear  problem 

(x'^x  +  yc'^Og  =  X^y,  (4) 

which  we  solve  numerically  using  Gaussian  elimination. 
The  computation  takes  less  than  20  seconds’  on  the  NBS 
mainframe  computer  for  N  =  128.  We  chose  7  =  0.01. 
The  solution  is  not  strongly  dependent  on  7  and  this 
value  produces  good  agreement  with  the  corresponding 
measured  waveforms,  where  these  are  available. 

We  test  the  algorithm  with  experimental  data  in 
this  manner.  Data  is  normalized  so  |Xj^|£1  and  |y|^|£1. 

Using  measured  data  x  and  y  ,  we  find  an 

s  s 

approximation  to  g,  g.  Typically  x^  is  a  step. 
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Having  thus  characterized  the  system,  we  input  a 
distinct  waveform  and  also  measure  the  output  y^. 

Using  g  and  y^,  we  calculate  an  approximation  to  x^, 

X  .  Typical  results  of  this  "double"  deconvolution 

are  discussed  below.  We  consider  the  agreement  to  be 
excellent.  The  selection  of  Y  is  not  addressed  here. 
In  a  future  paper,  we  will  present  a  detailed 
numerical  study  of  the  relation  between  Y  and  the 
level  of  noise. 

Experimental  Results 

We  apply  the  deconvolution  methods  detailed  above 
to  the  evaluation  of  several  sensors,  using  an 
experimental  configuration  which  has  been  described 
previously  [it].  Sensors  were  mounted  in  a  terminated, 
50-ohm,  coaxial  test  line  with  signal  input  from  a 
mer cury-wetted-relay ,  cable-type  pulser.  Signals 
propagate  through  the  line  past  the  sensor  location, 
resulting  in  the  application  of  different  transient 
waveforms  to  sensors  in  a  geometry  simulating  that 
encountered  in  pulse  power  or  gas- i nsul ated  high 
voltage  equipment.  A  variety  of  sensor  configurations 
can  be  accommodated  in  the  line,  which  is  constructed 
in  modular  form  to  allow  future  studies  of  current 
sensors.  For  the  work  described  here  two  sensors  were 
examined.  One  was  a  capacitive  sensor  with  a  coaxial 
low-side  capacitor  geometry.  The  second  was  a  model 
resistive  divider,  which  was  constructed  in  a  way  so 
that  its  response  characteristics  could  be  varied  by 
changing  the  effective  stray  capacitance.  As 
indicated  below,  it  was  operated  in  an  under¬ 
compensated  mode,  producing  an  overshoot.  Data  are 
acquired  with  a  transient  recorder  having  UOO-Mhz 
bandwidth  and  an  effective  digitizing  rate  above 
1  GHz. 

The  two  waveforms  used  in  the  present  studies  are 
shown  in  Fig.  1.  The  step-like  signal  is  used  in 
calculations  of  the  step  response  of  the  system,  while 
the  second  "test"  waveform  is  used  in  assessing  the 
quality  of  the  deconvolution.  The  test  waveform  was 
generated  by  inserting  a  small  inductive  element  in 
series  with  the  pulser.  These  waveforms  are  recorded 
directly  at  the  digitizer  and  represent  the  signals 
introduced  into  the  test  line.  These  direct  waveforms 
were  recorded  with  low  charge  voltages  on  the  pulser, 
which  resulted  in  an  erratic  behavior  for  the  test 
waveform  after  about  40  ns,  as  indicated  on  the 
figures  below.  At  the  normal  operating  voltages  used 
to  produce  signals  on  the  line,  this  was  not  observed. 
The  step  (and  impulse)  response  of  the  sensor  is 
calculated  using  eq.  4  with  the  direct  step  from  Fig  1 
and  the  output  from  the  sensor  as  input  and  output 
signals  respectively.  The  output  signals  from  a 
capacitive  sensor  for  the  two  input  waveforms  is  shown 
in  Fig.  2.  Using  the  calculated  step  response  and  the 
measured  output  for  the  test  waveform  input,  the  input 
signal  can  be  deconvolved,  again  using  eq.  4.  Fig.  3 
compares  the  input  signal  with  the  deconvolved  input 
signal,  where  the  two  curves  have  been  offset  to 
facilitate  comparison. 

A  similar  approach  is  followed  for  the  small 
resistive  divider.  The  undercompensated 
characteristics  of  this  sensor  are  indicated  in 
Fig.  4,  which  compares  the  divider  output  with  the 
test  waveform  input.  Deconvolved  and  input  waveforms 
for  this  divider  are  shown  in  Fig.  5. 

It  is  clear  that  there  is  good  agreement  between 
the  deconvolved  and  input  waveforms.  Serious 
deficiencies  in  the  response  of  a  practical  divider 
can  be  compensated  for  in  the  deconvolution  process. 
Similar  results  have  been  obtained  for  other  fast 


measurement  components,  including  an  E-dot  sensor  and 
a  long  signal  cable.  Methods  for  providing  a 
quantitative  method  for  assessing  the  quality  of  the 
deconvolution  are  being  investigated,  as  are  frequency 
domain  methods  like  those  discussed  in  [5]. 

Conclusions 


The  use  of  deconvolution  methods  in  evaluating 
sensors  and  in  compensating  for  imperfect  response  in 
practical  dividers  shows  considerable  promise. 
Additional  work  remains  to  be  done  to  assess  errors 
associated  with  the  method  before  full  implementation 
in  a  laboratory  measurement  system. 
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Figure  1 .  Waveforms  used  in  studies  measured  directly 
from  pulse  generator.  The  test  waveform 
does  not  replicate  beyond  40  ns  as 
indicated. 
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- DECONVOLVED  WAVEFORM 

- INPUT  WAVEFORM 


Figure  2.  Output  signals  from  a  capacitive  sensor  for 
test  line  input  signals  shown  in  Fig.  1. 


TIME  (ns) 

Figure  5.  Deconvolved  and  directly  measured  input 
signals  as  obtained  using  the  resistor 
divider.  Signals  have  been  offset  for 
comparison. 


Figure  3.  Comparison  of  deconvolved  and  input  signals 
for  the  capacitive  sensor.  Signals  have 
been  offset  for  comparison. 


Figure  4.  Comparison  of  the  output  signal  from  a 
resistor  divider  with  the  input  signal. 


